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ABSTRACT 

The variability of the high-energy gamma ray sources in the Third EGRET catalog is analyzed by a new 
method. We re-analyze the EGRET data to calculate a likelihood function for the flux of each source in each 
observation, both for detections and upper limits. These functions can be combined in a uniform manner with 
a simple model of the flux distribution to characterize the flux variation by a confidence interval for the relative 
standard deviation of the flux. The main result is a table of these values for almost all the cataloged sources. As 
expected, the identified pulsars are steady emitters and the blazars are mostly highly variable. The unidentified 
sources are heterogeneous, with greater variation at higher Galactic latitude. There is an indication that pulsar 
wind nebulae are associated with variable sources. There is a population of variable sources along the Galactic 
plane, concentrated in the inner spiral arms. 

Subject headings: methods: data analysis - methods: statistical - gamma rays: observations 



1. INTRODUCTION 

Of the 271 point sources of high-energy gamma-rays de- 
tected by the Energetic Gamma Ray Experiment Telescope 
(EGRET) on the Compton Gam ma Ray Observatory space- 
craft, most remain unidentified ( iHartman et alJll999l) . The 
nature of these sources is one of the biggest questions raised 
by the EGRET mission. Certainly some of them will turn 
out to belong to the established classes of sources: Active 
Galac tic Nuclei (AGN) ( Sowards-Emmerd et al. 2003 : | Willis! 
U996) and pulsars (Merck et al ..l996;.Mukheriee et alJll995t 
lYadiaaroglu & Romani 1995). A few, as well, may be discov- 
ered to be artifacts of the Galactic diffuse emission. However 
there are more unidentified sources than expected from these 
populations ( Merck et al. 1996). The large error regions for 
the EGRET sources often make it difficult to discern which 
of several candidates is the true counterpart in another wave- 
length. 

When studying the unidentified sources, there are sev- 
eral observables which one might use to determin e their 
nature . On e can loo k at t he spatial distribution ( Grenieri 
119951 iKan bach et al] 11995ft to estimate what fraction of 
the population is Galactic in nature, and what scale 
height the Galactic fraction has. One can look at po- 
sitional correla tions with other obj e ct types and Galac- 
tic structures ('Esposit o et"al] Il99 g: 'Gehrels et al.' '200Ct 
[Grenier 2000; Hartman et alJ199 ^ Kawasaki & Totani 2002; 
iMattox etalJ ll997; Reim er et alTl2003t iRomero et aljil99^ 
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iSUirner & Dermer 1995; Thompson et alJll99 ^ iTorres et alJ 
12003^ .Walker et aU .2002: ■Wallace'etan HoOB. The 
ergy spectrum is another characteristic which can 
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used to classify sources as "pulsar-like' 
"other" ( .Merck et al. 19961) . 

In this same vein, one can examine the EGRET catalog 
sources for evidence of time variability, in hopes of distin- 
guishing the various source classes. The known pulsars are 
seen to have a fairly constant flux (± < 10% when averaged 
over many pulse periods, consistent with the systematic un- 
certainties in EGRET flux measurements), as is expected from 
the nature of their energy production. Many AGN are seen 
to flare dramatically, which makes their identification easier. 
But, as with many aspects of the EGRET data, it is difficult 
to characterize the time variability of most sources, because 
of the limited statistics and the non-continuous observations. 
Thus it is natural to look for a statistic which will measure the 
variability in a rigorous way. 

This work is based on the dissertation of Tompkins (19991). 
Changes in the details of the method have caused changes in 
the numerical results, but the conclusions are not altered. 

In ^we briefly review the previous methods which have 
been used for evaluating the variability of EGRET point 
sources. In ^we outline an improved method. In ^we de- 
scribe how this method was applied to the EGRET data. The 
results are described in |5]and summarized briefly in |6] 

2. PREVIOUS METHODS 

All previous analyses of the variability of EGRET 
sourc es have used the flux values published in the Sec- 
ond (iThompson et all Il995h or Third EGRET Catalog 
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jHartmanetalJl999h . hereafter called 2EG and 3EG. This in- 
troduces the problem of upper limits. Many of the individual 
observations produced no positive detections. For these obser- 
vations, the flux is reported as a 95% confidence upper limit 
produced by LIKE, the standard EGRET likelihood program 
(iMattox et al. 1996). 

The first comprehensive method used for estimating the 
variability of EGRET sources is the one introduced by 
[McLaughlin et al. (1996). This method finds the of the 
measured source intensities in 2EG, assuming a constant flux, 
and uses V =-log2, where Q is the tail probability of obtain- 
ing such a large x^. There are two problems with this method: 
one statistical, and one conceptual. The statistical problem is 
introduced when using upper limits: McLaughlin et al. ( 1996) 
use zero as the flux estimate when the measured flux estimate 
would be negative, and for the error they use the LIKE upper 
limits. Thus their resulting "x^" will not have a x^ distribu- 
tion. Since the LIKE results are statistically conservative, this 
means that sources with upper limits included in the analysis 
will have a lower V than they should. The conceptual prob- 
lem with the V method is the use of a x^ statistic to determine 
variability. The Q statistic is the p-value for the hypothesis of 
constant source flux. Thus, sources with a large V are incon- 
sistent with being constant. But we do not know whether such 
a source has a large V because of large intensity fluctuations, 
or because of small error bars on the intensity measurements. 
Similarly, a source with a small V might truly be constant, or 
it might just have very poor measurements of its flux. 

Wallace et al. (2000) also applied this method to two-day 
segments of EGRET data, looking for rapid variation that 
would not have been discovered by examining full viewing 
periods. They were severly limited by large statistical uncer- 
tainties, but they found previously-undiscovered variation in 
s everal source s. 

iTorres et a n (I2001albft describe a statistic called the / in- 
dex. / is a normalized ratio of the standard deviation to the 
average flux. This is somewhat similar to the approach we de- 
scribe below. It is a step toward dealing with the conceptual 
problem described above, but it introduces another ad-hoc ap- 
proach to dealing with upper limits. Also there is no estimate 
of confidence intervals for / values. 

Other, less quantitative criteria have been used to produce 
lists of sources which are not highly variable. These lists 
are heterogeneous collections of sources with rather differ- 
ent properties. In addition to truly non-variable sources, they 
tend to include bright variable sources. They exclude mainly 
the sources which flared o nce or t wice and were otherwise 
undetectable. For instance. iGrenien ( 2000 ) defined a class of 
88 persistent unidentified sources. These are the ones which 
are significantly detected (> 4ct or > 5cr, depending on posi- 
tion) i n the cumulative data up to 1995 October Gehrels et al. 
(|2000) produced a similar list of 120 steady, unidentified 
sources. They required that the flux in the cumulative data 
must be either a more significant detection than any individual 
observation or within 3(T of the largest individual flux. Again 
the effects of average brightness and variability are hard to 
disentangle. It is difficult to compare sources or populations 
on the basis of these lists because they provide only a binary 
variable/nonvariable classification. 

3. A LIKELIHOOD APPROACH 
3.1. What Do We Want to Measure ? 



We propose to answer two questions about the variability 
of a source: How much does the flux vary and how precisely 
do we know this? To do this we define 

5^-. (1) 
M 

where /i and a are the average and standard deviation, respec- 
tively, of the true flux of the source. This fractional variabil- 
ity is much closer to what we really want to know about a 
source. There could also be other statistics which might be 
more meaningful for flaring sources, such as the peak flare 
flux divided by the quiescent level. Since the flux measure- 
ments have uncertainties, we must find a way to estimate 5 
from the data. Unfortunately this can't be done in a com- 
pletely model-independent way. 

3.2. Modeling the Source Flux Distribution 

We ignore the ordering and spacing of the individual flux 
measurements. We assume only that all flux values of a par- 
ticular source are random values drawn from the same proba- 
bility distribution. This approach elminates all consideration 
of long timescale correlations in the data. The probability 
distribution is parameterized by 5. Such a model is necessary 
in order to disentangle the effects of real flux variation from 
statistical fluctuations. An infinite number of source flux dis- 
tributions exist which will yield the same 5; we want to pick 
one that is compatible with our notions of what the true source 
flux distribution should be, but which is also fairly general. 

The first source flux distribution one might try is a Gaus- 
sian. Thus, we would find the likelihood of obtaining our data 
given that the true flux S was drawn from a Gaussian with 
mean /i and standard deviation a = fiS. We could then find the 
maximum likelihood value of 5, together with a confidence 
interval for 6 defined in the standard way (Eadie et al. 1971). 
Alternatively, we could follow the Bayesian procedure: form 
priors for the distributions of and 5, and find an estimate of 
5, together with an error region, by marginalizing over /i. The 
use of the Gaussian distribution for source fluxes has a flaw, 
however It allows the possibility of /i = 5 = oo), as well 
as the unphysical possibility of negative /i. 

Instead of a Gaussian, we choose to use the gamma distri- 
bution; that is, the flux S is drawn from the differential proba- 
bility distribution 

ca-l -S//3 

P(S)= , , (2) 

where the adjustable parameters a(> 0) and /3(> 0) describe 
the shape and scaling of the function, respectively. The pa- 
rameters can be related to the desired quantities by /i = /3a, 
a = f3^/a, and S = a"'/^. This distribution function has a peak 
if a > 1. For < a < 1, it decreases monotonically from 
S = 0. This behavior is a good match for the different types of 
observed variations. Some sources are steady (large a), while 
m any others are m ostly faint with occasional flares (small a). 

lTompkinsl ( II99 9) did the same analysis described below us- 
ing a lognormal distribution. There is an excellent correla- 
tion between the results, but the gamma distribution produced 
fewer numerical difficulties, such as lack of convergence. We 
take this as a validation of our choice of probability distribu- 
tion. 

3.3. Characterizing the Single-Measurement Likelihood 

It is standard in EGRET analysis to calculate t he likelihood 
of the observed data for a certain source flux (iMattox et alJ 
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Fig. 1 . — The likelihood curve for a sample source, together with a Gaus- 
sian of the same width, and a curve of the form in equation t3l . 



Il996h . Typically, the values output by an analysis program 
are the most likely flux and a confidence interval. In the case 
of upper limits, the most likely flux is taken as zero, and the 
confidence interv al might be defined in one of several ways 
jCasoetal.l1998h. 

For this work, we would like to be able to calculate the 
likelihood for any source flux, not necessarily one near the 
maximum likelihood value. Thus it is useful to look for a 
parameterized family of curves which will closely fit the full 
likelihood function. The results of the many individual mea- 
surements can be stored in a compact form and reproduced 
quickly as needed. 

After examining many such curves, we found a simple 
function which fulfills our requirements. We choose 



-(k+S/S,) 



(k+s/s,y 



(3) 



where 5,, k, and p are adjustable parameters and only S >Q 
is allowed. A likelihood function need not be normalized. If 
p > k, C{S) has a peak at S = {p- k)Sr of width S^^^/p, and 
the EGRET test statistic (Mat tox et aLl lT996) TS = 2p{k/p- 
ln(k/p)- 1). If p < k, C(S) has its maximum at 5 = 0. This 
function was originally derived by analogy with the Poisson 
probability function. The parameter p plays the role of an 
effective number of detected photons, while 1 /S'* is the expo- 
sure and the offset k accounts for the effects of diffuse emis- 
sion and nearby point sources. £(S) would become singular if 
^ < 0, but this never occurred in our analysis. 

This form fits the observed likelihood functions rather well 
(Figure[0, both in cases of large source flux and in upper limit 
situations. 

3.4. Calculating Likelihood 

With the parameterized distribution of source fluxes from 
the preceding section, we can use Bayes's Theorem to write 
the likelihood of observing a given viewing period's worth of 
data (Dj) from a source with a certain fi and S: 

/>oo 

£{fi,d,Dd= / dSiP{Di\Si)PiSi\fi,5) (4) 

JO 

dSiC(SdPiSi\fi,S) (5) 



dSi e 



-(ki^SilS.i) 



(Ai(52)l/5^r(l/52) 



^(l/5--l) -5,7(^52) 



(6) 



where the second term in equation (|6|l is the gamma distribu- 
tion written in terms of [i and 5. 

The likelihood of a sequence of observations (P\ ...Dm) 
will then be 

M 



or 



ln/:(M,'^) = f]ln/:(M,^,A). 



(7) 



/■=1 



3.5. Dealing with Time Scales 



We have assumed that the source flux during all the mea- 
surements was drawn independently from the same distribu- 
tion. This can be problematic if the measurements last for 
varying lengths of time. For example, in the EGRET data, 
there are frequently several successive viewing periods of a 
few days duration, pointing in nearly the same direction. This 
is in contrast to the normal mode of pointing in one direction 
for two weeks at a time. We do not want our results to depend 
on whether or not a two week viewing period was chopped 
into several bits. 

Our solution is to combine successive measurements if they 
occur within 30 days, and treat them as one measurement. In- 
stead of the single-observation likelihood in equation Q, we 
use a composite C(S) which is the product of the individual 
terms. The generalization of equation (|6j is obvious. These 
combined measurements are reflected in the number of obser- 
vations reported in Table 1 . 

Thus this analysis is sensitive to variations with a time scale 
of a month or mo re. We will no t be able to detect the 1-2 
day flares found bv lWallace et al J (120001) from several EGRET 
sources. 

4. ANALYSIS 
4.1. Method 

The first step in this analysis was the re-analysis of the 3EG 
catalog with a modified likelihood program which produces 
a likelihood function in the form of equation Q ( Tompkin^ 
Il999l) . For each Viewing Period, any of the 271 sources from 
3EG Catalog which were within the instrument's field of view 
were included in the source list. In addition, any of the 145 
marginal sources (which were used in the original analysis of 
the catalog but were not significant enough to merit inclusion 
in the official list) which were in the field of view were added 
as well. The maximum likelihood set of source fluxes was 
then found. For each source, the parameters characterizing 
£(5') described in ^3.3l were found. All source measurements 
were of flux (> 100 MeV), assuming a spectral index of -2.0, 
using all photons within 15° of each source, and including all 
nearby sources (both catalog and marginal) in the model. We 
use only data taken with the full EGRET field of view, which 
includes everything through Viewing Period 429.0 (ending 
1995 September 27) except for four periods (403.0, 403.5, 
411.1, and 41 1.5). With the smaller field of view used to take 
the rest of the data, we are not confident that the detector re- 
sponse functions are calibrated well enough for our purpose. 
Also the number of point sources detected in each observation 
is reduced, so the remaining data are less useful. Results are 
used only for the point sources within 25° of the center of the 
field of view. This data selection is nearly the same as that 
used to produce the 3EG Catalog. 
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The likelihood of the sequence of obser vations as a func- 
tion of S was then calculated, as in ^3.41 All observations 
within one month of each other were lumped together, as per 
^3.51 From this likelihood, the most likely values of S were 
obtained, as well as confidence intervals. The mean flux /i can 
be marginalized away as a nuisance parameter after checking 
that its most likely value is consistent with the 3EG flux. 

4.2. Sources Omitted from This Analysis 

In the results reported below, several 3EG sources are omit- 
ted for various reasons. 

Seven sources near the bright Crab and Vela pulsars 
are omitted because they are believed to be artifacts 
jThompson et al..2001J : 3EG J0824-4610, 3EG J0827-4247, 
3EG J0828-4954, 3EG J0841-4356, 3EG J0848-4429, 3EG 
J0859-4247, and 3EG J0521H-2147. 

The solar flare 3EG J0516H-2320 is omitted because it is 
clearly a unique occurrence. 

The sources 3EG J0829H-2413 and 3EG J1424H-3734 each 
had only one observation within 25° of the center of the field 
of view. Thus the variability of these sources can't be evalu- 
ated. 

A cross-check of this analysis and the 3EG flux values 
turned up three sources with serious discrepancies: 3EG 
J0556H-0409, 3EG J0628h- 1 847, and 3EG J 1 3 1 0-05 17. These 
three sources are omitted. The latter two are near very bright 
sources (Geminga and 3C279), so we believe that our modi- 
fied likelihood program may become confused in such areas. 
It is encouraging that only these three sources fail this test. 
The consistency of results for all the others encourages us that 
the modified likelihood program is valid in the other cases. 

4.3. Systematic Errors 

This analysis has two main sources of systematic er- 
ror. The first is from the instrument itself. The sensitivity 
of the EGRET i nstrument is variable (Hartman et al. 1999; 
iThompson et alJ>1995.) . Despite efforts to measure and com- 
pensate for lon g term drifts, there i s a residual variability of 
about 10-11% fesDo sito et alJl999l) . In an effort to minimize 
this systematic error, only source observations within 25° of 
the instrument axis were used in the analysis. 

It is conceivable that the variation in measured fluxes might 
be reduced by using the "noisy" measured sensitivity rather 
than the smoothed values that are usually employed in EGRET 
analysis. We tested this by calculating the fluxes of the 
three bright pulsars (Crab, Geminga, Vela) by both methods. 
With the measured sensitivity values. Vela appeared much 
more variable, the Crab slighly more, and Geminga slightly 
less. Thus it appears to be appropriate to continue using the 
smoothed sensitivities. 

Because of the systemitic variability, the Galactic and 
isotropic background levels Gm and Gb were allowed to vary 
from viewing period to viewing period for a given source in 
the likelihood analysis. Otherwise, a small change in a strong 
background could lead to a large change in a source flux. 
This proced ure is standard practice for EGRET data analysis 
dMattox etal.J996). 

The background levels were also allowed to float when CiS) 
was calculated. Unfortunately, computational limitations pre- 
cluded allowing nearby sources to vary during this calcula- 
tion as well. Thus, nearby sources were fixed at their max- 
imum likelihood intensities (calculated for that viewing pe- 
riod). This will lead to an underestimation of the error in flux. 



particularly for sources which have another source very close 
by. 

As both systematic effects will tend to exaggerate source 
variability, it is expected that even constant sources will show 
some small variability. To the extent that the systematics are 
dominated by the instrumental effects, the measured variabil- 
ity will be equal to the true variability added in quadrature to 
a constant systematic variability of about 10%. 

4.4. A Simple Variability Index 

Although the most likely value of 5 and its confidence inter- 
val provide a good description of a source's variability, there 
is considerable interest in the simpler question of whether a 
source is variable at all, independent of the question "how 
much?" This is the goal of McLaughlin et al. (1996). This 
problem is complicated by the presence of the systematic er- 
rors. To help answer the question, we define the statistic Vn- 
It represents the confidence with which we can reject the hy- 
pothesis (5 < 0. 12, where 0. 12 is a conservative estimate of the 
average systematic error 



yi2 = -logi„(l-P^2(r|l)), 



where 



r = 21n 



'Cmax(Mi'^ = 0-12)' 



(8) 



(9) 



£tnax(/^,<5) is the unconstrained maximum likelihood, 
Cms^(pL,5 = 0.12) is the maximum likelihood if 5 is con- 
strained to equal the systematic value 0.12, and P^2(r|l) is 
the cumulative probability distribution with one degree 
of freedom. The use of the d istribution can be justified 
by Wilks's Theorem dCashlfT979l) when several independent 
likelihood functions are combined, even though they may not 
be Gaussian functions individually. If the most likely 5 is less 
than 0. 12, then there is no estimate for Vn. Vn values greater 
than 16 are shown as oo in Table 1. The p-value for rejecting 
the hypothesis that a source is non-variable is 10"^'-. 

5. RESULTS 

The results of this analysis are shown in Table 1. The 
columns in the table give the names in the 3EG catalog, the 
most likely Galactic longitudes and latitudes, the most likely 
values of 5, the 68% lower and upper confidence limits on 5, 
Vn, the number of observations used in the analysis, possible 
associations with source classes, and specific identifications. 

5.1. Variability by Source Class 

As can be seen from Table 1, the limits on the variabil- 
ity of an individual source are usually not very strict. By 
grouping similar sources together, one can learn more about 
the average behavior of a class of sources. There are many 
sources which are found in multiple classes. Many of the 
unidentified sources are found in densely-populated regions 
containing various types of young-population objects. The 
error boxes of the EGRET sources often contain several plau- 
sible counterparts. Table 2 shows the results for the various 
source classes. It contains a description of each class, the 
number of sources in the class, the mean of 5, and an estimate 
of the intrinsic RMS dispersion of 5 in excess of the statisti- 
cal errors. The last two quantities are calculated by assuming 
a Gaussian distribution for 5 and finding the likelihood as a 
function of me an an d standard deviation by a method similar 
to that used in 
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Source Class Members (5) RMS(<5) Symbol 



3EG Classification 



Pulsars 


6 


0.11 ±0.02 


< 0.07 


PSR 


SNR associations 


10 


A 0"7 1 A AA 


< 0.19 


SNR 


ACjN ( A 111 3hCj) 


67 


A ir\ 1 A AO 


A n 1 A AC 

U.i/ ± 0.05 




Quasars 


C 1 

J 1 


n T? _L n no 

\j. 1 1 ± u.uy 


O.zj it O.Oo 




BL Lac objects 


15 


A yl A 1 A * 

0.49 ±0.16 


< 0.32 


BLL 


U maentmea, \b\ < 5 


47 


A /n 1 A Az; 

U.42 ibU.Uo 


AIT 1 A AO 

0.1 / ± O.Os 


UO 


Umaentined, 5 < \b\ < 15 


37 


r\ c 1 All 

0.56±0.12 


0.23 ± 0.09 


U5 


Umaentinea, 15 < o| < JO 
Unidentified, |/!)| > 30° 


45 


0.42 ±0.12 


< 0.16 


U15 


29 


0.84 ±0.15 


<0.27 


U30 


Possible isolated neutron stars 










Geminga-like objects 


3 


0.13±0.11 


<0.16 


Gem 


Pulsars without PWN 


7 


0.21 ±0.15 


<0.15 


noPWN 


Pulsars with PWN 


6 


0.66 ±0.29 


0.28±0.16 


PWN 


Gould's Belt, unidentified, \b\ > 5° 


33 


0.49±0.13 


<0.26 


G 


Persistent, unidentified (Gehrels et al.^OOff) 


88 


0.32 ±0.05 


0.10±0.05 


Per 


Steady, unidentified (Grenier 2000) 


120 


0.35 ±0.05 


0.11 ±0.04 


Std 


Associations from Romero et al. (1999) 










SNR 


17 


0.28 ±0.10 


<0.17 


RSN 


WR stars 


6 


0.45 ±0.22 


<0.27 


WR 


Of stars 


4 


0.26±0.15 


<0.19 


Of 


OB associations 


22 


0.40 ±0.10 


0.17±0.08 


OB 



To begin, consider four principal source classes based on 
the identifications in the 3EG catalog: Unidentified sources, 
Pulsars, Active Galactic Nuclei (AGN) (labeled "A" in 3EG), 
and sources which are spatially coincident with Supernova 
Remnants (SNR). The Unidentified source class was divided 
according to the Galactic latitude. These source classes 
were determined from the "source ID" (and by the "Other 
Name" category for the SNR associations) in the catalog. 
Also, 3EG J1048-5840 (PSR B1046-58) has been added as 
a sixth member of the pulsar class based on its recent detec- 
tion (|Kaspi et al. 2000). Altho ugh pulsations have bee n de- 
tected from PSR B1951H-32 ( Ramanamurthv et al. 1995), and 
reported from PSR B0656H-14 (Rama namurthv et al. 1996), 
they do not appear as detectable point sources in the 3EG cat- 
alog, so they are not considered here. 

As seen in Table 2, the pulsars, AGN, and unidentified 
sources clearly differ in their variability. The variability of 
the six pulsars is about 11%, not significantly higher than the 
predicted systematic variability for a constant source of 10%. 

When the firmly identified quasars and BL Lac objects are 
considered separately, a difference is clear. On the average, 
the BL Lacs are much less variable than the quasars. This 
confirms the observation made by Mukheriee et al. ( 1997). 

The ten sources in the 3EG SNR source class are some- 
what less variable than the low-latitude unidentified sources 
as a whole, and show significantly more variation than the 
pulsar class. Of course, it is by no means certain that these 
sources are supernova remnants. They are found in crowded 
areas of the sky along with ot her types of object s character- 
istic of star-fo rming regions ( Yadigar oglu & Romam 1997 ). 
iRomero et alJ 0999) identify 17 possible SNR associations, 
of which 7 are also in the 3EG SNR list. These two SNR col- 
lections have very similar variability properties. If the gamma 
rays are produced in a large region of the SNR, then the flux 
should be steady. Observable variation would be an indica- 
tion of gamma production in a pulsar wind or perhaps of the 
presence of a background blazar 

The change in unidentified source variability with latitude 
can be seen in Table 2. It is clear from the spatial dis- 



tribution of the unidentified sources that there are multiple 
types of astrophysical objects which have not been identi- 
fied (Gehrels et al. 2000; Grenier 1995; Kanbach et al. 1995). 
It is clear that the different types of unidentified sources 
have different variabilities as well. The high latitude sources 
have a variability index consistent with that of the identified 
AGN, and many of them have plausible AGN counterparts 
(iMattox et al. 2001, 1997; Sowards-Emmerd et al. 2003). The 
low latitude unidentified sources exhibit some variability, at a 
level a little more than the sources coincident with SNR. 

In the literature a number of other possible associations 
with source classes have been suggested, not based on the 
classifications in the 3EG catalog. Table 2 also summarizes 
t he proper ti es of these classes. 

iGrenied ll2003t) listed 13 3EG sources with positions 
consistent with radio pulsars and reasonable L^/E ra- 
tios. Six of these (3EG J1013-5915, 3EG J1420-6038, 
3EG J1837-0423, 3EG J1856H-0114, 3EG J2021H-3716, 
3EG J2227H-6122) have pulsar wind nebulae (PWN) which 
might be significant gamma sources. (3EG J1410-6147 is 
not counted as a PWN source because it can be associ- 
ated with two pulsars, but the one with a PWN can be dis- 
counted on energetic grounds.) The seven non-PWN sources 
(3EG J1014-5705, 3EG Jl 102-6103, 3EG J1410-6147, 
3EG J1639-4702, 3EG J1714-3857, 3EG J1824-1514, 
3EG J1837-0606) have an average 5 of 0.21 ±0.15, while 
the six PWN sources have an average 5 of 0.66 ± 0.29. These 
numbers suggest that most of the associations are correct and 
that PWNe can be strong, variable gamma emitters. Three 
of the PW N so urces are among the variable Galactic plane 
sources in ^5.41 

Four of the identified gamma-ray pulsars are known to 
have wind nebulae: Crab, Vela, PSR B 1046-58, and PSR 
B 1706-44. Their gamma emission is clearly dominated by 
the pulsars. Any variation in the nebular output in the > 100 
MeV band is hidden from this analysis. The Crab Nebula, fo r 
instance, is variable on a scale of years I de Jager e t al.'1996»). 
but the changes can be seen only in a certain energy band, 
and only in the off-pulse emission. Note: We do not con- 
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sider the complex source 3EG J0222+4253 here because of 
the c onfusion of emission from its pulsar and blazar counter- 
parts ( lHartmanetalJl999h . 

Three other 3EG sources have plausible associations 
with isolated neutron stars. They have hard spectra, 
large gamma to X-ray flux ratios, and no known pul- 
sar counterparts, s o they might be Ge minga-class objects: 
3EG JOOIOh-7309 ("Brazier et al.' 1998"), 3EG J1835H-5918 
jChandler et al. 2001; Mirabal & Halpern 2001; Reimer et al. 
lloOl), and 3EG J2020H-4017 ( Bra zier et al.l >1996). Individu- 
ally they show little sign of variability. As a class, their aver- 
age S is 0.13 ± 0. 11, quite consistent with the identified pul- 
s ars. 

lGrenied(|2000h s uggested an association of EGRET sources 
with Gould's Belt. iPerrot & Gr enier ( 2003) produced a model 
of the Belt based on a 3D simulation of Galactic dynamics. 
There are 47 unidentified 3EG sources whose positions are 
consistent with their model. Of th ese, four are omitted from 
this analysis as described in ^4.21 Ten more are at \b\ < 5°, 
so they are likely to be contaminated by the Galactic plane 
population(s). The remaining 33 have an average S of 0.49. 
This is a moderately low value, but these sources as a class are 
definitely more variable than the pulsars. Perhaps a fraction 
o f them are Myr - old pu lsars (|Hardin2 & Zhan2 2001). 

iRomero et alJ (1199 9^ made lists of possible associations 
with very young objects: WR stars. Of stars, and OB asso- 
ciations. These all have an average S of 0.26 to 0.45, fairly 
steady but more variable than pulsars . 

It has been pro posed (IColafrancescol l2002t 
iKawasaki & Totanil l2()02h that some high-latitude EGRET 
sources my be galaxy clusters which produce gamma rays in 
their intracluster gas. Of course such a s ource would be non- 
variab le on a scale of months to years. IKawasaki & Totanil 
( I2002h picked a sample of seven EGRET sources, cho- 
sen on the basis of steady emission and other criteria. 
3EG J13 10-05 17 is omitted from our analysis, and the other 
six show no significant evidence of variation. The upper 
limits on S are not tight (0.86 to 1.30), so the non-variability 
of these sources is an open question. Colafrancesco (2002) 
has a list of nine probable associations, three of which are 
also in the Kawasaki & Totani (2002) list. However, he 
includes 3EG J0253-0345 and 3EG J0215H-1123. Both of 
those are variable. 

As described in ^ Gehrels et al. (2000) produced a list of 
120 "steady" sources and iGrenien (i2000 ) a list of 88 "per- 
sistent" sources. Most of the "persistent" sources are also in 
the the "steady" list. These lists were attempts to eliminate 
strongly flaring sources and ones which are only marginally 
detectable. Apparently they succeeded, since the mean 6 val- 
ues (0.32 and 0.35) are fairly low. However, there are several 
examples of clearly variable sources in these lists. 

5.2. The Lowest Variability Sources 

With the variability of the source classes in hand, it is possi- 
ble to see which sources can be excluded from a source class 
based on the variability data. We begin by examining those 
sources which are inconsistent with being strongly variable: 
that is, the upper bound (Smax) of the 68% region for 6 is 
small. There are 8 sources for which this upper bound is less 
than 0.4, significantly lower than the typical value for AGN 
or unidentified sources. 

Five of the six pulsars are in the list, with PSR B 1055-52 
barely excluded. The possibility that others on the list might 
also be pulsars is intriguing. The sources are discussed indi- 



vidually below. 

The position of the Crab pulsar 3EG J0534H-2200, at 
the top of the list, is a bit surprising given its listing in 
McLau2hlin et al. (1996) as moderately variable. There are 
two reasons for this. T he first is the difference in the meth- 
ods: McLaughlin_eLalJreport the confidence level with which 
the hypothesis of steady emission can be rejected. In their 
analysis, the bright Crab, Geminga, and Vela pulsars are all 
inconsistent with being constant at the 95% level due to sys- 
tematic errors and the small statistical uncertainties. The sec- 
ond reason is the difference in data used: McLaughlin et al. 
used all measurements up to to 30° from the instrument axis. 
In this analysis, we have used only data out to 25°, because of 
the larger systematic errors at higher inclinations. 

3EG J1048-5840 is coincident with PSR B1046-58, a high 
E/d^ pulsar Although pulsed emission was not seen in the 
first three years of EGRET data (Fierro 1995), there has been 
a detection in more recent data iKasoi et alj200QI) . 

3EG J0222H-4253 is coincident wifli the BL Lac 3C 66A, 
and is 1° from the position of PSR J0218H-4232, from 
which pulsed emission above 10 MeV ha s been reporte d 
(Kuiper et aLl l200l iKuiper et aU EoOO; Verb unt et al.lll996h . 
In the 3EG catalog, the pulsar was not detected as a separate 
source, and it is suggested that the flux between 100 MeV and 
1 GeV is primarily from the pulsar, while the flux above 1 
GeV is mostly from 3C 66A (Hartman et al. 1999). The flux 
measurements of this "source" will be dominated b y the lower 
energy emission, hence the low variability. W allace et alJ 
(HbOO) found a significant variation of this "source" on a 2- 
day timescale, likely from the AGN. Since our analysis uses 
only full observations, we can't confirm or reject this varia- 
tion. 

The S value of the Vela pulsar 3EG J0834-451 1 is slightly 
higher than the 10% systematic variation. This is due to the 
presence of several artifact sources nearby (Thompso n et alJ 
2001). In a few viewing periods, the measured flux of Vela 
is low, and the flux from an artifact is quite high. Because 
nearby source fluxes were not re-fit for each point on the 
likelihood-vs-intensity curve, such "flux leaking" can lead to 
an apparently high variability. By fixing the fluxes of the spu- 
rious sources one could lower Vela's variability, but only if 
those nearby sources were not themselves variable. This was 
not done because of the unknown biases that could result. 

3EG J2020H-4017 is unidentified, but its position is coinci- 
dent with the 7-Cyg supernova remnant, so it is thus an in- 
triguing pulsar candidate. No pulsed signal has been detected 
from it, how ever, and the possibilit y exists that it is a radio- 
quiet pulsar dBrazier et al.il 996: .Chandler et al.l2001l) . 

3EG J 1734-3232, another unidentified source, is the last 
member of this group. 

It is noteworthy that 3EG J0721H-7120 almost qualifies for 
this list. It is ide ntified in the 3EG Ca t alog as the BL Lac 
objec t 07 1 6-H7 1 4 (iHartman et aP 1 19991 Ivon Montignv et alJ 
II995t rMukherie e et al.ll997h . The upper end of the 68% con- 
fidence region for S is 0.38, well below the typical AGN vari- 
ability of ^ 0.7. In large part, this is because this source had 
no dramatic flares during any observations. In addition, the 
viewing periods with the highest and lowest fluxes were each 
combined with less extr e me o bservations by the one-month 
averaging system ( ^4.11 ^3.5> . A more sophisticated model 
might pick up longer or shorter term trends. 

5.3. The Highest Variability Sources 
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Fig. 2. — Values of S, a measure of variability, for gamma sources within 
±5° of the Galactic Plane. For each source, the filled circle is the most likely 
value of 5 and the lower end of the bar is <5,^in, the 68% confidence lower 
limit. The vertical dashed lines are drawn at longitudes of ±55° to guide 
the eye. The horizontal lines are at the arbitrarily chosen cutoff values for 
variable sources: 0.7 for the most likely S and 0.2 for (5|^i„. All of the highly 
variable sources are within the central region of the Galaxy. 



From Table 1, 18 sources have a 68% confidence region 
with a lower bound greater than 0.7. These include 11 
sources identified as AGN. Four of the unidentifed sources 
(3EG J 1227+4302, 3EG J 1607-1 101, 3EG J2006-2321, 
and 3EG J2251-1341) are at high G alactic latitudes, and 
are gen erally believed t o be AGN jHartman et all Il999t 
iMattox et al. 2001; Walla ce et alJ2002h . The three remaining 
low latitude, unidentified sources are 3EG J 1704-4732, 3EG 
J1828+0142, and 3EG J1837-0423. Halpern et al. (2003) and 
[Spwards-Emmerd et al. ( 2003) propose a blazar identification 
for 3EG J1828+0142. .Tavani et aL C1997.) discoverd the vari- 
abihty of 3EG J 1837-0423, but they reject a blazar identifi- 
cation. None of the seven unidentified sources is classified as 
steady or persistent. 

Because the V statistic used by McLaughlin et al. ( 1996) in- 
dicates how improbable it is that the observations come from a 
constant source, the use of the 6 statistic in this way is nearly 
equivalent to the use of V to find variable sources. Indeed, 
when there is a close correspondence between 2EG and 3EG 
sources, sources with a high V also have a high lower limit on 
6. 

5.4. Galactic Plane Sources 

It is interesting to examine the behavior of the gamma 
sources near the Galactic plane. There are many which show 
no sign of variability, but there are some variable sources 
which are distributed non-uniformly in Galactic longitude. 
Figure 12] shows S values for the 66 sources within 6° of the 
Galactic plane. If we define a variable source as one with a 
most likely 6 value > 0.7 or a minimum S > 0.2, then it can 
be seen that almost all of the 17 variable, low-latitude sources 
are clumped within 55° of the Galactic center All 17 are 
unidentified. 

Figure|3]shows histograms of the most likely S values for all 
the sources within 6° of the Galactic plane. The dashed line 
shows the S distribution for the sources with \£\ < 55°, while 
the solid line shows the others. A Kolmogorov-Smirnov test 
dEadie et al . 1971) rejects the hypothesis that these two sets 
of values are drawn from the same population with a p- value 
of 5.1%. 
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Fig. 3. — Distributions of most likely S values for sources within 6° of the 
Galactic plane. Dashed line: Sources within 55° longitude of the Galactic 
center Solid line: Other sources. 



Figure|2]shows values of 5 for these sources as a function of 
longitude. The filled circles show the most likely S values and 
the ends of the bars show 68% confidence lower limits {Smin 
in Table 1). The horizontal dashed lines are drawn at S = 0.7 
and (5min = 0.2 to guide the eye. 

Some of the apparent variation might be due to flux leak- 
age between close companions, as with the Vela pulsar ( ^5.21 . 
Among the 17 variable sources there are three pairs and a 
triple within a few degrees. No pair of these sources has corre- 
lated fluxes with a correlation coefficient > 0.18, except 3EG 
J1823-1314 and 3EG J1826-1302. Their correlation coeffi- 
cient is -0.41, with a p-value of 2.4%. The variation of these 
two sources might be spurious, but the others seem to be on 
more solid ground. 

The distribution of variable sources is consistent with a 
young Galactic population, as trac ed by molecular cloud s 
(iDame et al."2001), OB associations fKaaret & Cottam> 1996h. 
SNRs and HII regions (Sturner & Dermer 1995; Sturne ret alJ 
1996; Yadiaaroalu & Romani 1997), WR stars and Of stars 
(Romero et al. 1999). The concentration in the inner radian 
and their latitude spread indicate typical distances of 5-8 kpc 
and luminosities L{E > 100 MeV, 47r sr) of (1-25) x 10^^^ 
(D /5 kpc)^ erg s~^ . This i s slightly higher than the range found 
by'Kanbach etaHfm?) for all the low-latitude sources in the 
2EG catalog, both variable and non-variable. 

This luminosity range is typical of X-ray binaries and 10 
kyr-old pulsars. The latter are expected to be steady emit- 
ters, but variable emission from pulsar wind nebulae is pos- 
sible. In particular, plerionic activity is more likely to ex- 
plain the hard but variable 3EG J1856+0114 source than 
pulsed emission from the 20 kyr-old pulsar or cosmic -ray 
emission in the surrounding W44 supernova shell. The 
same may be true for 3EG J 1420-6038 in the Kookaburra 
nebula and 3EG J1837-0423 near the G27. 8+0.6 plerion. 
3EG J 1809-2328 also coincides with a pulsar wind neb- 
ula, but inverse-Compton emission from electrons scatter- 
ing the nebular synchrotron emission or the ambient ra- 
diation field from the near by OB associatio n fail to ex- 
plain the source brightness ( Braie et all |2002^ ■ A 21 kyr- 
old pulsar, PSR J1826-1334, found near the edge of the 
3EG J 1826- 1302 error box, could power the gamma-ray 
source. Pulsars are also found near 3EG J1308-6112 and 
3EG J 1928+ 1733, but they are much too weak to explain the 
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Variability over months is common in accreting systems. 
One high-mass X-ray binary, GX 304-1, hosting a Be star 
and a neutron star on an eccentric 132.5-day orbit, is found 
toward 3EG J1308-61 12. It is known to exhibit strong long- 
term variability in X rays. Electrons accelerated at the shock 
between the pulsar and stellar winds have been proposed to 
shine up to TeV energies by upscattering the stellar radia- 
tion JKirk et al.1 1999; Tavani & Arons 1997), but the 272 s 
pulsar in GX 304-1 is not likely to power a strong, en- 
ergetic wind. 3EG J 1736-2908 coincides with the X-ray 
source GRS 1734-292, which has beenrecently identified as a 
nearby Seyfert 1 galaxy (Z = 0.0214) with radio jet emission 
and hard X-ray emission up to 400 keV (Marti et al. 1998). 
Confirmation of variable emission above 100 MeV from this 
galaxy would yield a first example of the activity of a Seyfert 
jet at high energy and suggest a valuable target for TeV tele- 
scopes. 

Curiously, 1 1 of the 17 variable sources are in the "steady" 
category, and 10 are "persistent." This happens because they 
are fairly bright. As long as they don't vanish entirely for 
many observations, they still meet the selection criteria. 

5.5. The Galactic Center 

The source nearest to the Galactic center, 3EG 
J1746-2851, seems to be one of these variable sources. 
iMaver-Hasselwande r et alJ ( Il998i) found only weak ev- 
idence that this source might be variable. With more 
data and a more sensitive method, the statistical evi- 
dence of its variability is fairly strong (V^ = 2.35). This 
would rule out any model in which this object is a cluster 
of neutron stars, a condensation o f dark matter parti - 
cles, or other non-compact object JBertone etal] I2OO3I 
Qilcaneo-Roldan & Moore 2000 
Markoff etalJll999t iMaver-Hasselwander 
1991 ! 

This conclusion is reinforced by Hoocer & Dingus ( 200^, 
who show that the position of the EGRET source is not con- 
sistent with the Galactic Center 3EG J 1746-2851 may be 
a field source like the other 16 low-latitude variables. The 
probablility of such a chance positional coincidence near the 
Galactic Center is about 2.7%. 

Since the diffuse emission is most intense in the neighbor- 
hood of the Galactic Center, it is conceivable that errors in the 
background model or instrument PSF will have their great- 
est effect there. To check this we correlated the flux of 3EG 
J 1746-2851 with the angle to the instrument pointing direc- 
tion. The correlation coefficient is 0.46, which has a p-value 
of 1.5%. If the maximum viewing angle is reduced from 25° 
to 17°, seven observations are eliminated. The flux-angle cor- 
relation is reduced to an insignificant level, but the estimate 
of 6 becomes consistent with zero. 

If there is a systematic problem at the Galactic Cen- 
ter, it should show up also in measurements of nearby 
sources. There are four other sources within 6° of the Center: 
3EG J1734-3232, 3EG J1736-2908, 3EG J1741-2312, and 
3EG J1744-3011. None of them shows such a large correla- 
tion of viewing angle and flux, and in one case the correlation 
is negative. 

However, 3EG J1741-2321 and 3EG J1744-3011 have a 
substantial negative flux correlation (-0.33 and -0.42 respec- 
tively) with 3EG J 1746-2851. This might be a sign of flux 
leakage, although 3EG J1741-2312 is rather far away. 

In view of these somewhat ambiguous signs of extra sys- 



1: iCesarinietalJ m 
;lwander et alJ ll998l1Poh]| 



tematic errors, we must declare that the variability of 3EG 
J 1746-2851 is only a tentative conclusion. 

6. CONCLUSIONS 

It is possible to characterize the flux variation of EGRET 
gamma-ray sources between viewing periods in a consistent 
way by using a simple model for the flux distribution and 
combining likelihood functions. This method produces use- 
ful results for almost all of the sources in the Third EGRET 
catalog, using both positive detections and upper limits in a 
uniform manner. The resulting table displays the average frac- 
tional variation 6 and its confidence interval. For faint sources 
the confidence interval is necessarily large, allowing a wide 
range of possible behaviors. 

The identified rotation-powered pulsars all show a small 
variability, consistent with the estimated 10% systematic er- 
rors in the flux measurements. In general, identified blazars 
are highly variable, although there are a few exceptions. The 
unidentified sources show a range of behaviors. Their av- 
erage variability increases with Galactic latitude, indicating 
that there are multiple populations. Unidentified sources as- 
sociated with pulsar with wind nebulae are, as a class, more 
variable than the ones associated with pulsars without nebu- 
lae. The highest-latitude unidentified sources have an average 
behavior similar to blazars. 

There is a population of 17 highly-variable sources along 
the Galactic plane. These are concentrated in longitude within 
1 radian of the Galactic center. This is probably a high- 
luminosity population of sources in the inner spiral arms and 
molecular ring, while the remaining less-variable sources are 
closer and less luminous. 
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Table 1 . Variability S of sources in the Third EGRET catalog 



3EG Name 


i 

fdeal 


b 

(deg) 


s 


S 

mm 


max 


Vi2 


Num 

Obs 


Classes 


ID 


JOO 10+7309 


119.92 


+10.54 


0.26 





0.73 


0.40 


3 


U5,SNR,Std,Pst,Gem 


CTA 1 


J0038-0949 


112.69 


-72.44 








1.34 




3 


U30,Std,Pst 




JO 11 8+0248 


136.23 


-59.36 


1.18 


0.48 


2.18 


1.64 


4 


a 


0119+041 


J0130-1758 


169.71 


-77.11 








0.63 




4 


a 


0130-171 


J0159-3603 


248.89 


-73.04 








0.85 




3 


U30,Std,Pst 




J0204+1458 


147.95 


-44.32 


1.29 


0.60 


2.54 


2.62 


3 


QSO 


0202+149 


J0210-5055 


276.10 


-61.78 


0.31 


0.15 


0.57 


1.26 


6 


QSO 


0208-512 


J0215+1123 


153.75 


-46.37 


1.27 


0.59 


2.48 


2.10 


3 


U30 




J0222+4253 


140.14 


-16.77 








0.37 




4 


BLL 


3C 66A, PSRJ0218+4232 


J0229+6151 


134.20 


+1.15 








0.51 




5 


U0,Of,OB,G,Std,Pst 




J0237+1635 


156.77 


-39.11 


0.89 


0.46 


1.89 


5.07 


3 


BLL 


0235+164 


J0239+2815 


150.21 


-28.80 








0.41 




5 


a 


0234+285 


J0241+6103 


135.87 


+0.99 


0.38 


0.17 


0.79 


1.41 


5 


U{),G,Std,Pst 




J0245+1758 


157.62 


-37.11 


1.14 


0.44 


2.16 


1.38 


4 


U30,Std 




J0253-0345 


179.70 


-52.56 


1.34 


0.63 


2.98 


2.03 


2 


U30 




J0323+5122 


145.64 


-4.67 


0.92 


0.31 


1.68 


1.26 


7 


U0,G,Std 




J0329+2149 


164.90 


-27.88 








1.55 




5 


U15,Std 




J0340-0201 


188.00 


-42.45 








1.24 




2 


QSO 


0336-019 


J0348+3510 


159.06 


-15.01 


0.45 





1.29 


0.24 


6 


U15,G,Std,Pst 




J0348-5708 


269.35 


-46.79 


1.19 





1.86 


0.79 


8 


U30 




J0404+0700 


184.00 


-32.15 


0.39 





1.53 


0.11 


6 


U30,Std,Pst 




J0407+1710 


175.63 


-25.06 


0.90 





1.59 


0.55 


9 


U15,G,Std 




J0412-1853 


213.90 


-43.29 








2.10 




2 


QSO 


0414-189 


J0416+3650 


162.22 


-9.97 


0.59 


0.10 


1.15 


0.89 


9 


a 


0415+379 


J0422-0102 


194.88 


-33.12 








1.08 




4 


QSO 


0420-014 


J0423+1707 


178.48 


-22.14 


0.42 





1.02 


0.42 


10 


U15,G,Std,Pst 




J0426+1333 


181.98 


-23.82 








0.73 




12 


U15,G,Std,Pst 




J0429+0337 


191.44 


-29.08 








0.73 




8 


U15,Std,Pst 




J0433+2908 


170.48 


-12.58 


0.40 


0.10 


0.75 


0.82 


11 


QSO 


0430+2859 


J0435+6137 


146.50 


+9.50 








0.54 




6 


U5,Std,Pst 




J0439+1105 


186.14 


-22.87 








0.74 




12 


U15,G,Std,Pst 




J0439+1555 


181.98 


-19.98 


1.08 


0.48 


1.63 


1.20 


12 


U15,G 




J0442-0033 


197.20 


-28.46 


1.59 


1.12 


2.31 


14.59 


8 


QSO 


0440-003 


J0450+1105 


187.86 


-20.62 


1.13 


0.78 


1.64 


8.42 


12 


QSO 


0446+112 


J0456-2338 


223.96 


-34.98 


0.21 





1.42 


0.05 


5 


QSO 


0454-234 


J0458-4635 


252.40 


-38.40 


0.41 





1.19 


0.23 


6 


QSO 


0454-463 


J0459+0544 


193.99 


-21.66 


0.74 





1.43 


0.54 


10 


QSO 


0459+060 


J0459+3352 


170.25 


-5.77 


0.59 


0.14 


1.18 


0.95 


12 


U5,G,Std,Pst 




J0500+2529 


177.18 


-10.28 








1.22 




12 


U5,G 




J0500-0159 


201.35 


-25.47 


1.09 


0.58 


1.71 


2.29 


9 


QSO 


0458-020 


J05 10+5545 


153.99 


+9.42 








0.41 




7 


U5,Std,Pst 




J0512-6150 


271.25 


-35.28 








0.71 




8 


a 


0506-612 


J0520+2556 


179.65 


-6.40 








0.53 




12 


U5,G,Pst 




J0530+1323 


191.37 


-11.01 


0.74 


0.55 


1.05 


oo 


13 


QSO 


0528+134 


J0530-3626 


240.94 


-31.29 


0.58 





1.75 


0.56 


3 


a 


0521-365 


J053 1-2940 


233.44 


-29.31 


0.85 





1.91 


0.60 


4 


a 


0537-286 


J0533+4751 


162.61 


+7.95 








0.51 




8 


U5,Std,Pst 




J0533-6916 


279.73 


-32.09 


0.29 





0.79 


0.13 


8 




LMC 


J0534+2200 


184.56 


-5.78 


0.08 


0.03 


0.15 




12 


PSR 


Crab 


J0540-4402 


250.08 


-31.09 


0.75 


0.46 


1.27 


6.18 


6 


BLL 


0537-441 


J0542+2610 


182.02 


-1.99 


0.57 





1.17 


0.76 


12 


U0,SNR,RSN 


S147 


J0542-0655 


211.28 


-18.52 


1.38 


0.84 


2.14 


1.93 


7 


a 


0539-057 


J0546+3948 


170.75 


+5.74 


0.10 





0.68 




11 


U5,Std,Pst 




J0613+4201 


171.32 


+11.40 


0.72 


0.09 


1.38 


0.87 


10 


U5,Std,Pst 




J06 16-0720 


215.58 


-11.06 








1.32 




8 


U5,G,Std 




J06 16-33 10 


240.35 


-21.24 


0.17 





1.04 


0.04 


3 


U15,G,Std,Pst 




J0617+2238 


189.00 


+3.05 


0.25 


0.10 


0.46 


0.82 


12 


UO,SNR,RSN,Std,Pst 


IC443 


J0622-1139 


220.16 


-11.69 


1.04 


0.26 


1.96 


1.08 


5 


QSO 


0616-116 


J063 1+0642 


204.71 


-1.30 


0.01 





1.16 




13 


U0,RSN,Std 




J0633+1751 


195.13 


+4.27 


0.10 


0.06 


0.16 




10 


PSR 


Geminga 


J0634+0521 


206.18 


-1.41 


0.01 





1.12 




13 


U0,Of,OB,Std 




J0702-6212 


272.65 


-22.56 


1.01 


0.34 


1.65 


1.06 


9 


U15,Std 




J0706-3837 


249.57 


-13.76 








1.75 




4 


U5,G 




J0721+7120 


143.98 


+28.02 








0.41 




5 


BLL 


0716+714 


J0724-4713 


259.00 


-14.38 


1.16 


0.57 


1.94 


2.02 


6 


U5,G,Std 




J0725-5140 


263.29 


-16.02 


0.89 


0.29 


1.64 


1.18 


7 


U15,G,Std,Pst 




J0737+1721 


201.85 


+18.07 








0.70 




3 


BLL 


0735+178 


J0743+5447 


162.99 


+29.19 


1.21 


0.67 


2.06 


4.78 


5 


QSO 


0738+5451 


J0747-3412 


249.35 


-4.48 


0.59 





1.55 


0.29 


6 


UO,WR,G,Std 




J0808+4844 


170.46 


+32.48 








0.65 




5 


a 


0804+499?, 0809+483? 


J0808+5114 


167.51 


+32.66 








1.18 




4 


a 


0803+5126 


J0808-5344 


268.24 


-11.20 


1.04 


0.55 


1.68 


2.16 


8 


U5,G,Std 




JOS 12-0646 


228.64 


+ 14.62 








0.90 




2 


a 


0805-077 


J0821-5814 


273.10 


-12.04 


1.16 


0.46 


1.84 


1.05 


8 


U5,G 




J0828+0508 


219.60 


+23.82 








2.17 




2 


BLL 


0829+046 


J0834-4511 


263.55 


-2.79 


0.16 


0.10 


0.28 


0.61 


7 


PSR 


Vela 
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3EG Name 


I 

(deg) 


b 

(deg) 


5 


<5 

min 


max 


Vl2 


Num 

Obs 


Classes 


ID 




J0845+7049 


143.54 


+34.43 


0.62 


0.23 


1.43 


1.38 


4 


QSO 


0836+710 




J0852-1216 


239.06 


+19.99 


1.21 


0.52 


2.26 


1.94 


4 


QSO 


PMN J0850- 


-1213 


J()853+1941 


206.81 


+35.82 








0.81 




3 


BLL 


0851+202 




J09()3-3531 


259.40 


+7.40 


0.52 





1.28 


0.70 


7 


U5,Std 






J0910+6556 


148.30 


+38.56 


0.64 





1.61 


0.60 


4 


U30,Std 






J0917+4427 


176.11 


+44.19 








0.53 




5 


a 


0917+449 




J0952+5501 


159.55 


+47.33 


0.39 





0.88 


0.31 


9 


QSO 


0954+556 




J0958+6533 


145.75 


+43.13 


0.96 


0.42 


1.72 


2.29 


6 


BLL 


0954+658 




J1009+4855 


166.87 


+51.99 








0.94 




9 


a 


1011+496 




J1013-5915 


283.93 


-2.34 


0.14 





0.52 


0.18 


9 


UO,G,RSN,Std,Pst,PWN 






J1014-5705 


282.80 


-0.51 








0.55 




10 


UO,G,Std,Pst,noPWN 






J1027-5817 


284.94 


-0.52 


0.21 





0.46 


0.39 


9 


U0,OB,G,Std,Pst 






11045-7630 


295.66 


-15.45 








0.44 




7 


U15,Std,Pst 






J 1048-5 840 


287.53 


+0.47 








0.31 




8 


PSR,Std 


PSR B 1046- 


-58 


J1052+5718 


149.47 


+53.27 


0.24 





1.02 


0.10 


6 


a 


1055+567 




J1058-5234 


285.98 


+6.65 








0.47 




8 


PSR 


PSR B 1055- 


-52 


J 1102-6103 


290.12 


-0.92 








1.08 




7 


U0,WR,OB,G,SNR,Std,Pst,noPWN 


MSH 1 1-62 




Jl 104+3809 


179.83 


+65.03 


0.35 





0.80 


0.46 


7 


BLL 


Mrk421 




Jl 133+0033 


264.52 


+57.48 


0.82 





1.51 


0.65 


8 


U30 






Jl 134-1530 


277.04 


+43.48 


1.12 


0.61 


1.84 


3.79 


7 


a 


1 127-145 




J 1200+2847 


199.42 


+78.38 


1.17 


0.75 


1.77 


5.49 


9 


QSO 


1156+295 




J1212+2304 


235.57 


+80.32 








1.13 




11 


U30 






J1219-1520 


291.56 


+46.82 


1.10 


0.56 


1.72 


1.77 


9 


U30 






J1222+2315 


241.87 


+82.39 


1.09 


0.60 


1.65 


1.61 


10 


U30 






J1222+2841 


197.27 


+83.52 


0.62 


0.22 


1.20 


1.17 


8 


BLL 


1219+285 




J1224+2118 


255.07 


+81.66 


0.45 





0.98 


0.62 


9 


QSO 


1222+216 




J1227+4302 


138.63 


+73.33 


1.36 


0.76 


2.37 


2.34 


4 


U30 






J1229+0210 


289.95 


+64.36 


0.46 


0.22 


0.85 


1.56 


8 


QSO 


30 273 




J1230-0247 


292.58 


+59.66 


0.59 





1.31 


0.71 


8 


QSO 


1229-021 




J1234-1318 


296.43 


+49.34 


0.61 





1.33 


0.77 


8 


U30,Std,Pst 






J1235+0233 


293.28 


+65.13 


0.21 





0.86 


0.11 


8 


U30,Std,Pst 






J 1236+0457 


292.59 


+67.52 


0.53 





1.42 


0.10 


8 


a 


1237+0459 




J 1246-0651 


300.96 


+55.99 


0.60 


0.11 


1.23 


0.89 


8 


QSO 


1243-072 




J1249-8330 


302.86 


-20.63 


0.71 





1.64 


0.62 


6 


U15,Std 






J1255-0549 


305.09 


+57.06 


0.90 


0.65 


1.35 


oo 


8 


QSO 


30 279 




J 1300-4406 


304.60 


+18.74 


0.54 





1.40 


0.30 


8 


U15,Std 






J 1308+8744 


122.74 


+29.38 








0.94 




6 


U15,Std 






J 1308-61 12 


305.01 


+1.59 


0.72 


0.29 


1.45 


1.37 


8 


U0,OB,G,Std 






J 13 14-3431 


308.21 


+28.12 








0.42 




6 


a 


1313-333 




J13 16-5244 


306.85 


+9.93 


0.42 





0.98 


0.49 


8 


U5,G,Std,Pst 






J 1323+2200 


359.33 


+81.15 


1.09 


0.52 


1.79 


1.61 


7 


a 


1324+224 




J1324-4314 


309.32 


+19.21 








0.49 




8 




Cen A 




J1329+1708 


346.29 


+76.68 


0.62 





1.36 


0.14 


9 


QSO 


1331+170 




J1329-4602 


309.83 


+16.32 








0.76 




8 


U15,Std,Pst 






J1337+5029 


105.40 


+65.04 


0.53 





1.30 


0.41 


5 


U30,Std,Pst 






J1339-1419 


320.07 


+46.95 


0.81 





1.56 


0.34 


7 


QSO 


1334-127 




J1347+2932 


47.32 


+77.50 


0.64 





1.60 


0.40 


6 


U30,Std 






J1409-0745 


333.88 


+50.28 


1.42 


1.00 


2.07 


oo 


8 


QSO 


1406-076 




J1410-6147 


312.18 


-0.35 


0.32 


0.08 


0.64 


0.76 


8 


U0,Of,OB,G,SNR,RSN,Std,Pst,noPWN 


G312.4-0.4 




J 1420-6038 


313.63 


+0.37 


1.03 


0.37 


1.83 


1.59 


8 


U0,OB,G,Std,Pst,PWN 






J 1429-42 17 


321.45 


+17.27 


0.93 


0.44 


1.58 


2.15 


9 


QSO 


1424-418 




J 1447-3936 


326.12 


+17.96 


0.28 





1.07 


0.09 


9 


U15,G,Std,Pst 






J1457-1903 


339.88 


+34.60 


0.49 





1.65 


0.07 


6 


U30,Std 






J 1500-3509 


330.91 


+20.45 








0.92 




9 


U15,G,Std,Pst 






J1504-1537 


344.04 


+36.38 


1.24 


0.58 


2.15 


1.79 


5 


a 


1504-166 




J15 12-0849 


351.49 


+40.37 








0.52 




5 


QSO 


1510-089 




J1517-2538 


339.76 


+26.60 








1.59 




6 


a 


1514-241 




J1527-2358 


342.97 


+26.50 


0.92 





1.71 


0.07 


9 


U15 






J 1600-0351 


6.30 


+34.81 


0.01 





1.62 




4 


U30 






J1605+1553 


29.18 


+43.84 








0.70 




2 


BLL 


1604+159 




J1607-1101 


0.91 


+29.05 


1.45 


0.93 


2.10 


0.93 


10 


U15 






J1608+1055 


23.03 


+40.79 


1.09 


0.25 


2.72 


1.15 


2 


QSO 


1606+106 




J1612-2618 


349.76 


+17.64 


1.04 


0.48 


1.62 


1.47 


12 


U15,G 






J1614+3424 


55.15 


+46.38 








0.41 




2 


QSO 


1611+343 




J1616-2221 


353.00 


+20.03 








0.54 




12 


U15,G,Std,Pst 






J1621+8203 


115.53 


+31.77 








0.52 




4 


U30,Std 






J1625-2955 


348.82 


+13.32 


1.62 


1.22 


2.19 


oo 


12 


QSO 


1622-297 




J 1626-25 19 


352.14 


+16.32 


1.05 


0.54 


1.60 


2.03 


13 


QSO 


1622-253 




J 1627-24 19 


353.36 


+16.71 








0.57 




12 


U15,G,Std,Pst 






J1631-1018 


5.55 


+24.94 


0.35 





1.01 


0.21 


12 


U15,G,Std,Pst 






J 163 1-4033 


341.61 


+5.24 


0.31 





0.90 


0.16 


14 


U5,G,Std 






J1633-3216 


348.10 


+10.48 


0.59 





1.25 


0.50 


13 


U5,G,Std 






J1634-1434 


2.33 


+21.78 








0.88 




12 


U15,G,Std,Pst 






J1635+3813 


61.09 


+42.34 








0.55 




2 


QSO 


1633+382 




J1635-1751 


359.72 


+19.56 


0.83 





1.48 


0.16 


13 


U15,G 






J1638-2749 


352.25 


+12.59 


0.47 


0.17 


0.88 


1.06 


14 


U5,G,Std,Pst 







12 



Table 1 . Variability 5 of sources in the Third EGRET catalog — Continued 



3EG Name 


i 


b 

(deg) 


s 


5 

mm 


max 




Num 

Obs 


Classes 


ID 


J1638-5155 


334.05 


-3.34 


0.19 





0.91 


0.13 


13 


UO,Std,Pst 




J1639^702 


337.75 


-0.15 








0.57 




12 


U0,OB,RSN,Stcl,Pst,noPWN 




J1646-0704 


10.85 


+23.69 


0.75 





1.41 


0.76 


12 


U15,G,Std,Pst 




J1649-1611 


3.35 


+17.80 








1.21 




14 


U15,G,Std,Pst 




J 1652-0223 


15.99 


+25.05 








0.96 




6 


U15,G,Std,Pst 




J1653-2133 


359.49 


+13.81 


0.95 





1.51 


0.73 


14 


U5,G 




J1655-4554 


340.48 


-1.61 


0.78 


0.18 


1.50 


0.97 


13 


U0,WR,OB,Stcl,Pst 




J1659-6251 


327.32 


-12.47 


0.73 





1.66 


0.47 


6 


U5,Std 




J1704-^732 


340.10 


-3.79 


1.37 


0.91 


1.93 


2.67 


13 


UO 




J1709-0828 


12.86 


+18.25 


0.67 





1.40 


0.37 


11 


U15,G,Std,Pst 




J1710-4439 


343.10 


-2.69 


0.07 





0.21 




14 


PSR 


PSRB1706^ 


J1714-3857 


348.04 


-0.09 


0.17 





0.52 


0.29 


13 


UO,RSN,Std,Pst,noPWN 




J1717-2737 


357.67 


+5.95 


0.91 


0.50 


1.43 


2.72 


15 


U5,Std,Pst 




J1718-3313 


353.20 


+2.56 


0.82 


0.43 


1.36 


2.50 


15 


UO.OB 




J1719-()430 


17.80 


+18.17 








0.81 




8 


U15,G,Std,Pst 




J1720-7820 


314.56 


-22.17 


0.64 





2.13 


0.02 


3 


a 


1716-771 


J 1726-0807 


15.52 


+14.77 


0.45 





1.03 


0.36 


12 


U5,G,Std,Pst 




J 1727+0429 


27.27 


+20.62 








0.74 




3 


QSO 


1725+044 


J1733+6017 


89.12 


+32.94 


0.29 





1.32 


0.15 


4 


U30,Std 




J1733-1313 


12.03 


+10.81 


0.37 


0.14 


0.68 


1.01 


12 


QSO 


1730-130 


J1734-3232 


355.64 


+0.15 








0.37 




16 


U0,OB,RSN,Std,Pst 




J1735-1500 


10.73 


+9.22 


0.87 





1.45 


0.38 


14 


U5,G 




J 1736-2908 


358.79 


+1.56 


0.64 


0.31 


1.12 


1.67 


15 


U{),Std,Pst 




J1738+5203 


79.37 


+32.05 


0.51 


0.10 


1.37 


0.85 


3 


QSO 


1739+522 


J1741-2050 


6.44 


+5.00 








0.57 




14 


U5,Std,Pst 




J1741-2312 


4.42 


+3.76 


0.26 





0.78 


0.15 


14 


U0,Std,Pst 




J1744-0310 


22.19 


+13.42 


0.59 





1.49 


0.52 


7 


QSO 


1741-038 


J1744-3011 


358.85 


-0.52 


0.30 





0.61 


0.64 


15 


U0,RSN,Std,Pst 




J1744-3934 


350.81 


-5.38 


0.60 


0.04 


1.14 


0.82 


12 


U5,Std,Pst 




J1746-1001 


16.34 


+9.64 


0.42 





0.88 


0.54 


12 


U5,G,Std,Pst 




11746-2851" 


0.11 


-0.04 


0.48 


0.29 


0.75 


2.35 


15 


U0,RSN,Std,Pst 




J 1757-07 11 


20.30 


+8.47 


0.51 





1.11 


0.42 


14 


U5,G,Std,Pst 




J1800-0146 


25.49 


+10.39 








0.58 




8 


U5,G,Std,Pst 




J1800-2338 


6.25 


-0.18 


0.16 





0.58 


0.26 


16 


UO,SNR,RSN,Std,Pst 


W28 


Jl 800-3955 


352.45 


-8.43 








1.13 




13 


QSO 


1759-396 


Jl 806-5005 


343.29 


-13.76 


0.88 





1.56 


0.34 


10 


a 


PMNJ1808-5011 


Jl 809-2328 


7.47 


-1.99 


0.71 


0.46 


1.13 


3.93 


15 


U0,OB,Std,Pst 




J1810-1032 


18.81 


+4.23 








0.58 




14 


UO.Std 




J1812-1316 


16.70 


+2.39 


0.82 


0.47 


1.35 


2.73 


15 


UO,Std,Pst 




J1813-6419 


330.04 


-20.32 








0.97 




4 


U15,Std,Pst 




J1822+1641 


44.84 


+13.84 


1.11 


0.54 


2.10 


2.26 


4 


U5,G 




J1823-1314 


17.94 


+0.14 


0.60 


0.21 


1.18 


1.18 


14 


U0,OB,SNR,Pst 


Kes67 


J1824+3441 


62.49 


+20.14 


0.01 





1.30 




8 


U15,G 




J1824-1514 


16.37 


-1.16 








0.84 




14 


U0,OB,RSN,Std,Pst,noPWN 




J1825+2854 


56.79 


+18.03 








1.33 




7 


U15,G 




J1825-7926 


314.56 


-25.44 


0.76 





1.77 


0.78 


4 


U15,Std,Pst 




J1826-1302 


18.47 


-0.44 


0.88 


0.50 


1.45 


3.22 


14 


U0,OB,Std,Pst 




Jl 828+0142 


31.90 


+5.78 


1.60 


1.02 


2.48 


6.89 


6 


U5 




Jl 832-21 10 


12.17 


-5.71 


0.62 


0.19 


1.20 


1.05 


13 


QSO 


1830-210 


J 1834-2803 


5.92 


-8.97 








0.59 




14 


U5,Std,Pst 




J1835+5918 


88.74 


+25.07 


0.15 





0.47 


0.09 


4 


U15,Std,Pst,Gem 




Jl 836-4933 


345.93 


-18.26 








0.97 




9 


U15,Std,Pst 




J 1837-0423 


27.44 


+1.06 


1.50 


0.97 


2.26 


2.95 


8 


UO,RSN,PWN 




J 1837-0606 


25.86 


+0.40 








0.67 




8 


U{),Std,Pst,noPWN 




J 1847-32 19 


3.21 


-13.37 


0.75 





1.36 


0.71 


14 


U5,Std 




Jl 850+5903 


88.92 


+23.18 








1.58 




5 


U15 




Jl 850-2652 


8.73 


-11.76 


0.81 





1.42 


0.59 


13 


U5 




J1856+0114 


34.60 


-0.54 


0.71 


0.28 


1.53 


1.57 


6 


UO,SNR,RSN,Std,Pst,PWN 


W44 


J1858-2137 


14.21 


-11.15 








0.80 




12 


U5,Std,Pst 




J1903+0550 


39.52 


-0.05 








0.42 




5 


UO,SNR,RSN,Std,Pst 


G40.5-0.5 


J1904-1124 


24.22 


-8.12 


0.08 





0.76 




10 


U5,Std,Pst 




J1911-2000 


16.87 


-13.22 


0.30 





0.82 


0.25 


11 


QSO 


1908-201 


J1921-2015 


17.81 


-15.60 


1.24 


0.80 


1.76 


2.19 


12 


a 


1920-211 


J1928+1733 


52.71 


+0.07 


0.76 


0.33 


1.57 


1.83 


6 


UO 




J1935^022 


358.65 


-25.23 


1.34 


0.83 


2.12 


4.41 


6 


QSO 


1933^00 


J1937-1529 


23.95 


-17.12 


1.35 


0.72 


2.05 


1.08 


8 


QSO 


1936-155 


J1940-0121 


37.41 


-11.62 


1.28 


0.70 


2.10 


2.48 


6 


U5 




J 1949-3456 


5.25 


-26.29 


1.17 


0.64 


1.84 


2.59 


8 


U15 




J1955-1414 


27.01 


-20.56 


0.84 


0.10 


1.62 


0.88 


7 


U15,Std 




J1958+2909 


66.23 


-0.16 


0.32 





0.86 


0.31 


9 


UO,Std,Pst 




J1958-4443 


354.85 


-30.13 








1.64 




4 


U30 




J1959+6342 


96.61 


+17.10 








0.76 




6 


U15,Std,Pst 




J2006-2321 


18.82 


-26.26 


1.31 


0.73 


2.12 


2.63 


6 


U15 




J2016+3657 


74.76 


+0.98 


0.44 





1.02 


0.63 


9 


U0,WR,OB,RSN,Std,Pst 




J2020+4017 


78.05 


+2.08 


0.06 





0.24 




9 


U0,OB,G,SNR,RSN,Std,Pst,Gem 


7Cyg 


J2020-1545 


28.09 


-26.62 








1.30 




4 


U15,Std 



Variability of EGRET Gamma-Ray Sources 



Table 1 . Variability 5 of sources in the Third EGRET catalog — Continued 



3EG Name 


i 


b 


s 


5 

min 


max 


Vi2 


Num 


Classes 


ID 






("deal 










Obs 






J2021+3716 


75.58 


+0.33 


0.36 


0.03 


0.83 


0.71 


9 


U0,WR,OB,Std,Pst,PWN 




J2022+4317 


80.63 


+3.62 


0.16 





0.80 


0.12 


9 


U{),WR,OB,G,Std 




J2025-0744 


36.90 


-24.38 


0.87 


0.42 


1.81 


3.14 


4 


QSO 


2022-077 


J2027+3429 


74.08 


-2.36 








0.43 




9 


U0,OB,Std,Pst 




J2033+4118 


80.27 


+0.73 


0.21 





0.52 


0.40 


9 


U0,Of,OB,G,Std,Pst 




J2034-3110 


12.25 


-34.64 


1.18 


0.55 


2.05 


1.43 


5 


U30 




J2035+4441 


83.17 


+2.50 


0.05 





0.66 




9 


UO,G,Std,Pst 




J2036+1132 


56.12 


-17.18 








1.30 




5 


BLL 


2032+107 


J2046+0933 


55.75 


-20.23 








0.80 




6 


U15,Std 




J2055-4716 


352.56 


-40.20 


0.93 


0.34 


2.09 


1.59 


3 


QSO 


2052-474 


J2100+6012 


97.76 


+9.16 


0.13 





0.70 


0.03 


6 


a 


2105+598 


J2158-3023 


17.73 


-52.25 


0.59 





1.66 


0.75 


3 


BLL 


2155-304 


J2202+4217 


92.59 


-10.44 


0.80 


0.02 


1.57 


0.85 


8 


BLL 


BL Lac 


J2206+6602 


107.23 


+8.34 








0.49 




5 


a 


2206+650 


J2209+2401 


81.83 


-25.65 


0.91 





1.86 


0.36 


5 


QSO 


2209+236 


J22 19-7941 


310.64 


-35.06 








0.89 




4 


U30,Std,Pst 




J2227+6122 


106.53 


+3.18 


0.20 





0.82 


0.21 


5 


U0,OB,G,Std,Pst,PWN 




J2232+1147 


77.44 


-38.58 


0.48 


0.18 


0.97 


1.22 


7 


QSO 


2230+114 


J224 1-6736 


319.81 


-45.02 








1.54 




3 


U30 




J2243+1509 


82.69 


-37.49 


1.06 





1.79 


0.86 


7 


U30 




J2248+1745 


86.00 


-36.17 


0.65 





1.48 


0.54 


7 


U30,Std,Pst 




J2251-1341 


52.48 


-58.91 


1.30 


0.73 


2.19 


4.04 


5 


U30 




J2254+1601 


86.11 


-38.18 


0.52 


0.31 


0.91 


3.93 


7 


QSO 


2251+158 


J2255+1943 


89.03 


-35.43 


1.18 


0.25 


1.95 


0.95 


7 


a 


2250+1926 


J2255-5012 


338.75 


-58.12 


0.50 





2.00 


0.33 


2 


U30 




J23 14+4426 


105.32 


-15.10 


1.11 


0.20 


2.39 


0.93 


3 


U15,Std 




J2321-0328 


76.82 


-58.07 


1.36 


0.78 


2.27 


2.46 


5 


QSO 


2320-035 


J2352+3752 


110.26 


-23.54 


1.37 


0.72 


2.43 


2.23 


4 


a 


2346+385 


J2358+4604 


113.39 


-15.82 








0.61 




3 


QSO 


2351+456 


J2359+2041 


107.01 


^0.58 


0.74 





1.57 


0.53 


7 


QSO 


2356+196 



Note. — The sources marked "a" in the Classes column are low-confidence AGN identifications in the 3EG catalog. For the other Classes, see Table 2. 



^ See text for comments about the variability of 3EG J 1746-2851. 



